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We give a detailed proof of the conjecture that the asymptotic entanglement cost of preparing a 
state p is equal to lim n ^oo Ef(p® n )/n where Ef is the entanglement of formation. 
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I. INTRODUCTION 



One of the central issues in quantum entanglement theory is to determine how to optimally convert between different 
entangled states shared by distant observers Alice and Bob. More precisely, one is interested in the conversion of 
m pairs of particles, with each pair in a state p, into n pairs, each in another state p' , by means of local quantum 
operations and classical communication (LOCC) [Q, so that the asymptotic ratio m/n is minimal. Of course, the 
perfect transformation by LOCC 

is usually impossible. Therefore, one permits imperfections and requires only asymptotically perfect transformations: 
the state p® m is transformed into some state p' n , which for large n approaches pf® n with a chosen distance measure 
D: 



lim D(p' n ,p'® n )=0 



If the final state p is a two-qubit singlet state ^ (|01) — 1 10)) (which we will denote by l^-)), then the process of 
conversion is called distillation If, instead, it is the initial state that is in the form of singlets, then one refers 
£^ ■ to formation In this paper we are interested in the latter process. We will call the optimal asymptotic yield the 
entanglement cost, and denote it by E c . In Ref. [|| it was shown that for a pure state p = \ip)(ip\, E c is equal to the 
entropy of either of its reductions, e.g. pa = Trs(p). Thus, to produce |<0)®™ one needs m w nS(pA) singlets, i.e. the 
initial state has to be \^_)® nS ^ A \ 

That result suggested the following stochastic method for the production of p out of singlets |Q . One decomposes p 
into an ensemble p — p^ipi) (ipi\ of pure states. Then one picks a state \ipi) according to the probability distribution 
{pi}, makes the state \ipi) from initially shared singlets, and finally forgets the identity of the state. Therefore, one 
needs on average piS(pf) singlets (where pf is the reduction of \tpi))- One can choose the ensemble which minimizes 
the above average, with the corresponding minimal cost being called the entanglement of formation of p JlJ , which we 
will denote by Ef(p). 

The above scenario can be improved if we realize that it might be more economical to produce the state p® n all 
at once than it is to produce its n constituents one by one (see e.g. Q). Thus, the proposed optimal cost should be 
revised to (5) 

E?(p)= lim -M—L. (1) 

(In the following we will argue that the limit on the right-hand-side of this equation exists.) This quantity is believed 
to be equal to the entanglement cost of preparing the state /?, E c (p). The definition of the entanglement cost, however, 
refers to a deterministic number of input singlets, while in the stochastic method the number of input singlets is a 
random variable. Furthermore it is not clear that the stochastic method is really the optimal way to produce the state 
p. In this paper we resolve these issues by proving that Ef°(p) is equal to the entanglement cost E c (p). Our result is 
in a sense dual to that of Rains |(| concerning the entanglement of distillation. He showed that the entanglement of 
distillation, if defined in a way analogous to our definition of the entanglement cost, is equal to the expected yield of 
a stochastic protocol for generating output singlets. 
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To begin, let us sketch our approach. In our proof, we will first show that Ey(p) > E c (p), i.e. that there exists a 
formation protocol that achieves the asymptotic rate ^ ~ Ey (p) , by explicitly constructing such a protocol based 
on the law of large numbers in both its classical and quantum forms. Next, we will show that Ey(p) < E c (p). 
The latter inequality will be derived using only some general properties of the entanglement of formation, such as 
its monotonicity under LOCC operations and its quite strong continuity properties, proved by Nielsen 0. Such an 
approach, focussing on abstract features of quantum information-theoretic quantities, has already proven to be fruitful 
in the domain of quantifying entanglement as well as in the study of quantum channel capacities [^],[l0) , providing a 
simple view of intrincate topics. To illustrate its power we will now sketch the proof of the inequality Ey(p) < E c (p). 

As mentioned before, the entanglement of formation is monotonically decreasing under LOCC operations Q], i.e. 

E f {p)>E f {C{p)), (2) 

for any state p and LOCC operation C. Moreover, the entanglement of formation is continuous Q in the sense that 
for states p and p' we have 

\E f (p)-E f (p')\ <5D(p,p')log 2 dimH + 2r,(D(p,p')), (3) 

where p and p' are supported on the Hubert space Ji and r\(x) = —xlog 2 x, under the assumption that D(p,p') is 
sufficiently small. D is the Bures distance, given by D(p,p') = 2yjl - F{p,p') with F(p,p') = Tr^p 1 / 2 p'p 1 / 2 (Note 
that the normalization factor for D is not completely standardized. We have made our choice to agree with Ref. ||). 
The function F is called the Uhlmann fidelity or square-root-fidelity |fl],|f2|. Inequality (||) implies that if two states 
are close to one another, then so are their densities of entanglement. A similar continuity result has been proved in 
Ref. [jl3| for the relative entropy of entanglement |l4| ]. 

Now, in order to prove Ey < E c , consider the optimal formation protocol, i.e. the sequence of LOCC operations 

A n producing states p n « p® n out of |^>_)® m , so that m/n tends to E c . If we take large n, so that D(p® n ,p n ) is 
small, then by continuity (|^) we have 

ZAP®") ^ Ef(Pn) (4) 

n n 

Since p n = A„(|'I , _)® m ), then from the monotonicity of Ef we obtain 



< " ~ ' = — (5) 



n 



Now, by defnition, the left-hand-side of the equation §Q) tends to EJ 3 while the last term of estimate (||) tends to E c , 
hence we obtained the required inequality. 

Since we used only two properties of the entanglement of formation, we can rephrase the result in a more general 

setting. Consider any function /, which can be regularized, i.e. for which the limit f°°(p) = lim.,^00 — - exists. 
Now, if / is monotonic and continuous in the sense of equation (^), then /°° is a lower bound for E c . This supports 
the view, according to which the dual measures of entanglement, entanglement of distillation Ed and entanglement 
cost E c are in a sense extreme ones || . 

II. ENTANGLEMENT COST AND ENTANGLEMENT OF FORMATION 

Let us now pass to the fully rigorous part of the paper. Analogous to the definition of asymptotic reducibility in 
Ref. jl5) (cf. Jl6|]), we define the asymptotic entanglement cost for the preparation of a bipartite state to be 



EAp) =inf{£MVe>0,<5>0, 3m,n,£, \E - — | < <5 and L>(£(|*_)(*_p m ), p®") < e} , (6) 

n 



where l^-) is the singlet state in C2 <S> C2, C is an LOCC superoperator and D is again the Bures distance. 
Our main result is the following: 

Theorem 1 The asymptotic entanglement cost of preparing a state p is given by 

EM . lim ace, p, 
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where E c (p) is defined in Eq. (Q) and Ef(p) is the entanglement of formation of p, defined as 

E f (p) = mm £>£(|ifcWi|), (8) 
£={pi,\it>i)} j 

and E(\ip)(tp\) is the von Neumann entropy of the reduced density matrix of 

To make sense of the above claim, we note that the sequence ^— ^ — has a well-defined limit. This is a 

consequence of the fact that if a sequence (a„) satisfies a n < cn for some constant c and a n + a m > a n+m for any 
m, n, then (a n /n) is convergent fll?|j . It is easy to see that our sequence satisfies the conditions. 

We begin by proving that E c obeys a form of additivity that is highly desirable for an asymptotic cost function, 
and that will be useful in what follows. 

Lemma 1 E c {p® k ) = kE c (p) for all k = 1, 2, . . .. 

Proof One direction is simple: since a protocol to approximate p® k is just a protocol to approximate k copies of p, 
E c (p® k ) > kE c (p). The first step in demonstrating that the inequality holds in the opposite direction will be to 
show that in the definition of E c , n can be taken to be arbitrarily large. Suppose not, in other words, that for fixed 
e,<S > 0, it is impossible to choose n> N,m, and £ such that \E c (p) - m/n\ < 5 and D(£(|*_)(*_ \® m ), p® n ) < e. 
By the definition of E c (p) and the fact that the set of LOCC operations is closed, this implies, however, that there 
exist n<N 7 m 7 and £ such that E c (p) = ni/n and £(|*_)(*_ \® m ) = p® n . Setting r = \N/n\ and applying £® r to 
|\E'_)( X I'_ |® mr then violates the assumption that n could not be taken arbitrarily large. 

Now, to complete the proof, we will pick n ^> k in order to efficiently approximate some number of copies of p® fc . 
Formally, we write n = rk + s, where r and s are non-negative integers and s < k. If we produce an approximation 
to within e of p® n starting from m singlets and then throw away the extra s copies of p, we are left with an 
approximation to within e of p® rk that still required m singlets. The waste, however, becomes insignicant since 
\E — m/n\ — * \E — m/rk\ as n — > oo. □ 

Note that full additivity of E c (p), i.e. E c (p®a) — E c (p) + E c (a) would be a stronger statement, possibly requiring 
additivity of Ef. With that technical lemma concluded, Theorem |l| becomes a consequence of the following two 
lemmas. 



Lemma 2 Let the entanglement of formation of a (finite dimensional) density matrix p be Ef(p) = 
where the optimal ensemble is given by {pi, \ipi)} k =1 . We have Ef(p) > E c (p). 

Proof In the limit of large n, we will approximate p® n by 

»i 

where jj 7 ^ is the (strongly) typical set defined as follows. All states \ tp s ) (unnormalized) in this set are tensor products 
of states y/pi\ipi) and they are such that every state ^Jpl^i) occurs pin ± Sin log p . (2)/fc times in the tensor product. 
(We assume that pi < 1 for all i. Otherwise, p is a pure state.) This implies that the total probability p s — {ip s \il) a ) 
for a state \ip s ) is bounded as 

2 -n(H(p))+8 1 ) < ps < 2 - n (H(?))-8l) : (10) 

where H{p) is the Shannon entropy of p = [pi). The density matrix p„(n) is basically the state p® n constructed 

by cutting off the unlikely sequences of states and then renormalizing |is|| . In order to bound the fidelity for this 
approximation, we add these unlikely sequences of states with probability in the sum. For any 5i > 0, t\ > and n 
sufficiently large we have by the asymptotic equipartition theorem |l!| that the total probability p T (n) for the typical 

set is larger than 1 — ei, where e\ — > as n — > oo. Thus, using the joint concavity of F pTfl, we can bound 



F(p® n ,p T ^)>VT^. (ii) 

Consider the density matrix p T (n) and its decomposition in terms of the states \ips}- Each state \tp s ) can be obtained 

from a set of EPR pairs by entanglement dilution H. In particular, let \ip s ) be a state in which every state \ipi) occurs 
Pin ± Sin\og p . {2)/k times. Starting from a set of maximally entangled EPR pairs, we do entanglement dilution 
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for each state |^ lglp *™ ± ' 5inloE Pi( 2 )/ fc . For any 62 and €2 greater than zero there exists an n such that starting from 
(pin + Sinlog pi (2)/k)[E f (\il}i}{ipi\) + S 2 ] EPR pairs we can obtain an approximation to |^.}®Pin±*i«l°s Pi ( 2 )/ fc which 
has square-root-fidelity larger than 1 — 62- Since there are k states in the optimal ensemble (and k is finite), we can 
therefore approximate the state \ij} s ) with square-root-fidelity > (1 — e 2) fe , starting from n J2iPiEf(\ipi}(if)i\) + 

n(0(6i) + 0(S 2 )) EPR pairs, with e 2 0, Si -> 0, and 5 2 for n -> 00. 

The approximation of p T <n) by p' (n) = ^ er (ro l^s) (V's 1 5 where |^) is the approximation of \ip s ) which we obtain 

by entanglement dilution starting from the set of EPR pairs, has the property that 

F( PtM ,p' tW )> £ K^W)I>(1- e 2 ) k = l-e 3 , (12) 

where £3 — > for n — > 00, since k is finite. Furthermore, since we can make every state \ip' s ) starting from a given 
set of EPR pairs, we can make any convex combination, for example p' ( „ } , of the states \ip' s ) (see [|o| and also pl|), 

starting from this same set of EPR pairs. 

Finally, we can use the triangle inequality for the Bures metric, and Eqs. ([ll]) and ( |l2] ) to obtain that 

D(p*»,(/ w) < 2^1-VT^T 1 + 2^T 3 . (13) 

which is the desired result. □ 

The lemma can be applied to bound E c (p) from above by Ef(p® k )/k where k is any fixed number of copies of p, 
using Lemma [l| Consequently, we have that 

Ef{p)>E c (p). (14) 
Let us now prove the converse of this relation. 
Lemma 3 Ef{p) < E e (p). 

Proof The basic idea of the proof, as sketched in the introduction, is to use the continuity of the entanglement 
of formation [|| and its monotonicity under LOCC operations. Throughout the proof we will use the notation 

A n = ZliZ > . Now suppose that the lemma does not hold, so that limA„ > E c (p). It follows that there exists an 
integer N, such that for all k > N, 

\A k -lunA n \<A = limA «- E ^ >0. (15) 

Let pk = jC ( I SI/ ){ 1 J / |® m ) be an approximation to p® k . Now, supposing that p is supported on Hilbert space H, fix e 

such that 5elogdim7i — 2eloge = A. It then follows from Eqs. (||) and our choice of e that if D(p® k , pt) < e then 

\E f (p® k )-E f ( Pk )\<kA. (16) 

Now we will apply the definition of E c {p) to fix k. From the definition of E Cl Eq. (||), and the proof of Lemma [j], we 
have that there exists a k > N as well as m and C such that 

\E c (p)-^\<A (17) 



and D{p® k ,p k ) < e, where p k = £(|*_)(*_ \® m ). 
Next, we can bound 



hW„-^ 



< |limA„ - A k \ 



< 2A, (18) 



by using Eqs. dH) and ©. 

This gives Ef(pk) > fc(lirnA n — 2A). On the other hand, by Eq. (fTi]), we have 



^ / (|*_X*-| 8ro ) = m < k(E c (p) + A), (19) 

which finally yields 

E f (p k ) - E f (\*_){*_\® m ) > k(]imA n - E c (p) - 3A) = kA, (20) 
a contradiction since the entanglement of formation cannot increase under LOCC operations. □ 
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III. ALTERNATIVE DEFINITIONS OF THE ASYMPTOTIC COST 



While E c (p) is perhaps the most natural function to associate with the asymptotic cost of preparing a bipartite 
state, other definitions would have been consistent with our discussion in the introduction. An example of a different 
but perhaps useful definition of the asymptotic entanglement cost is the following: 



E alt (p) = itf{E\Ve>0,6>0, 3N \ (Vn > N 3m, £, \E-— I < (5 and £>(£(!*_) I®" 1 ), p® n ) < e)\ . (21) 

n 



The difference between E c and E a u is that for fixed fidelity and cost the former only requires the existence of a 
single n such that p® n can be approximated efficiently while the second requires that there exist a threshold N such 
that for all n > N the state p® ra can be approximated efficiently. One consequence of this difference is that it is 
immediately clear that E a i t is additive in this sense of Lemma [l], while some work was required to prove that E c was. 
Indeed, a priori, one might not expect the two definitions to agree. Suppose, for example, that for fixed e there exists 
a protocol to make an approximation p n of p® n from m singlets, with the approximation satisfying D(p n , p® n ) = e. 
In order to produce an approximation to p® nk from km singlets, one might think that applying this protocol k times 
would be a good strategy. One finds, however, that for e > 0, 



lim D(p® nk ,p® k ) = lim 2 Jl - (1 - e 2 ) fe = 2. (22) 

k — >oo k — *oo V 

The protocol, therefore, generally fails to produce a good approximation for large k. This example suggests that the 
existence of a threshold N beyond which approximations to some given fidelity are always possible is a very difficult 
condition to satisfy. Nonetheless, by applying the results of the previous section, it is actually easy to see that the 
definitions E c and E a it are equivalent, so that the extra condition can be met without increasing the asymptotic unit 
cost. First, notice that the argument of Lemma || actually also works for E a it so that E a it(p) < E^ip)- Next, since 
the definition of E a i t is more stringent than that of E c , we have that E c {p) < E a i t (p). Combining these inequalities 
with the result of Lemma |3[ we get 

Ef{p) < E c (p) < E alt (p) < Ef(p), (23) 

so that these two definitions of the entanglement cost always agree. 

A final pair of alternative definitions for the asymptotic entanglement cost would use the trace distance 

d( P) <7) = ^Tr\p-a\, (24) 

in place of the Bures distance. Lemma |l| is easily seen to hold for d since its proof uses only axiomatic properties of 
all metrics and stability with respect to tensor products but no other properties specific to D. The status of Lemmas 
U and |3| is slightly more involved. The trace distance and Bures distance are related by the following inequalities |2^] 



-D(p, of = \- F(p, a) < d(p, a) < ^l-F(p,ay = sj — - — , (25) 

which show that d and D are equivalent metrics and bounds them in terms of each other by functions that are 
independent of the dimension of the underlying Hilbert space, 7i. In other words, d and D are equivalent metrics even 
in the asymptotic regime. To be concrete, applying the right-hand inequality to Eq. (|l3|), for example, is sufficient 
to show that Lemma || also holds for the trace distance. In order to prove Lemma pL we used the continuity relation, 
Eq. (||), but effectively only required the weaker inequality 

\E(p)-E(a)\ < BD(p,a) logdimft + C*, (26) 

where B and C are constants and Ti is the supporting Hilbert space. Applying Eq. (|25|), however, gives an inequality 
of the form 



\E(p)-E(a)\ <ByJd{p,a) logdimW + C, (27) 

which again is sufficient to carry through the rest of the proof. Therefore, all our conclusions hold even if the 
entanglement cost is defined using the trace distance. Indeed, they should hold for any metric equivalent to the Bures 
metric in which the equivalence is given by a function that is independent of the dimension of TL. 
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IV. CONCLUSIONS 



We have given two rigorous definitions of the asymptotic cost of preparing a bipartite mixed state p and shown 
them both to be equal to the regularized entanglement of formation, £JS°(p), resolving an important conjecture in 
the theory of quantum entanglement. Furthermore, we have shown that this asymptotic cost is fairly insensitive to 
the choice of metric on density operators. In particular, the Bures distance and trace distance result in identical 
asymptotic costs. 

An important problem left open by this work is the question of actually evaluating Ey (p) . Even the non-regularized 
function Ef(p) is notoriously difficult to calculate; its value is only known for some very special cases |Q,^3|- If it turns 
out that Ef is not additive for tensor products, then, in spite of the results of this paper, determining the asymptotic 
cost of preparing a state remains quite a formidable problem. 
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